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1. Introduction 

In the last decade basically the complete automation of QCD Nlo calculations has been 
achieved [1]. The predictions of these calculations depend on the QCD input parameters, such 
as the renormalisation scale /^r, the factorisation scale /Tp, the strong coupling CCs{lJ.^) and the 
parton density functions (PDFs) fa{x,lJL^) for the gluon and the (anti-)quarks. 

These input parameters need to be varied when we want to calculate theory uncertainties, e.g. 
when we follow the convention to vary the scales with factors 0.5 and 2 to estimate the effect of 
higher-order corrections. But for QCD cross section calculations beyond the Lo it usually takes 
days to achieve a precision that is comparable to data of current experiments, e.g. at the Large 
Hadron Collider. So a full re-evaluation for a given set of input parameter values uses up lots 
of CPU cycles. This problem is even more severe for another application of parameter variation: 
Fits of ttv and PDFs, where many re-evaluations are needed before the fit converges. Certainly for 
PDFs, which have many degrees of freedom in a fit, performing dedicated calculations for each 
PDF candidate value becomes impossible. 

The solution is to identify how a prediction depends on the QCD input pai'ameters. Then we 
can reuse the independent and CPU-intense parts for subsequent re-evaluations. This has been 
achieved for those types of QCD calculations, where the input parameters (or a known function of 
them) enter the calculation as prefactors to the identified independent parts [ 2 ]. 

There are two approaches to re-evaluate these independent parts for another choice of input 
parameter values. One is direct reweighting: We present here a novel feature of the 2.2.0 release 
of the Monte-Carlo event generator Sherpa [3] for performing an on-the-lly reweighting. Another 
approach is to create an interpolation grid, which allows for even faster re-evaluations. Currently 
there are two implementations for those grids, APPLgrid [4] and FAStNLO [5], and two interfaces 
for their automated creation using NLO events, aMCfast [6] and MCGRID [7]. Here we discuss 
the new MCGRID 2.0 release and present its newly added support for FAStNLO and for filling the 
fixed-order expansion of Mc@Nlo calculations. We also present a method for quantifying the 
dependency of a parton shower algorithm on the QCD input parameters. This is a first step towards 
the goal of correctly taking care of these dependencies in a reweighting or an interpolation grid 
creation. 


2. QCD NLO Events 


Consider the structure of an Nlo QCD calculation using Catani-Seymour subtraction [8], 
which consists of the following building blocks: The Born term (B), the virtual correction plus 
I-operator of the integrated subtraction terms (VI), the K- and P-operators of the integrated sub¬ 
traction terms including the collinear factorisation remainders (KP), and the real-emission (R) and 
correlated subtraction (D 5 ) terms: 
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with A/triai = Lf ^triai i the sum of the number of trials for each event and with the observable cuts 
and projections 0{^), where <I> gives the phase space point. The 0{^) are only dependent on the 
final state and so are independent of the QCD input parameters. We will now make explicit the 
dependency structure wich respect to the input parameters of each of the building blocks B, VI, KP 
and D 5 . The simplest are the “Born-like” terms X = B, R, D 5 , where a combination of QCD input 
parameters occurs as a total prefactor, as is well-known from the Lo factorisation theorem: 

X{^)=faixa,^j)Mxb,^j) (<!>), ( 2 . 2 ) 


where n is the number of strong vertices in the LO diagrams and p is either 0 (B) or 1 (RjD^). 
Here and in the following the primed quantity (X') on the right has no dependency on the QCD 
input parameters. We call it the weight of the corresponding term (X). As one can see, for the 
evaluation of the PDFs the longitudinal momentum fractions Xa/b have to be known in addition to 
the weight. The structure of the VI term is more complicated due to the occurring renormalisation 
scale logarithms: 


VI(<I>b) =fa{Xa,lll) fb{Xb,l4) 


<+' (At«) VT(<I>b) + {cfh + \ 


(2.3) 


with coefficients and scale logarithms Ir = log where }Xr is the scale at which the had 
been evaluated in the first place. So for re-evaluating the VI term, not only the weight VT and the 
momentum fractions x^jb need to be provided, but also the coefficients and the original scale 
flR. The last missing building block are the KP terms: 


+ fa{Xa,I^F) 


(2.4) 


with being the ratio of x^/b before and after the initial-state branching, and the coefficients 
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for / = 0,..., 3 with 


If = log^. 
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(2.5) 


The Cp and Cp are independent of the QCD input parameters. Again, jlp is the scale at which the 
cf had been evaluated. The functions give the correct combination of PDFs for the possible 
initial-state branchings, depending on whether a quark or a gluon is split, cf. [7]. 

All of these weights, coefficients and kinematics information can be passed in a standard 
HepMC [9] event record: The record comes with a HepMC : : WeightContainer object, which 
can be used to store all weights and coefficients. This must of course be supported by the event 
generator, as done by e.g. Sherpa [3]. 

With this information provided, the weights and coefficients can be combined with any arbi¬ 
trary choice of QCD input parameters. The same is true for the Nlo expansion of Mc@Nlo and 
(partly) for multi-leg merging calculations, but their structure will not be discussed here to keep 
this note short. It is not true however for parton showers, their inherent all-orders resummation 
structure renders such a direct decomposition to prefactors composed of input parameters and a 
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fixed number of weights/coefficients impossible. This lack of a correct treatment for variations in 
the shower results in a caveat for PDF variations of multi-leg merging calculations: They are not 
fully exact, because for them the PDFs enter via the shower branching probabilities determining 
the rejection weights for emissions already at the matrix element level. 


3. Sherpa-Internal Reweighting 

The newly added internal reweighting feature of Sherpa 2.2.0 uses the extraction of weights 
and coefficients sketched in the previous section to perform an on-the-fly reweighting for vari¬ 
ations of the QCD input parameters. Any number of combinations of scales and PDFs can be 
specified prior to the run. The results are inserted into the HepMC : : WeightContainer object 
and thus can be either saved to disk, programmatically accessed or directly passed via the internal 
interface to the Rivet analysis framework [10]. Lo, Nlo, the fixed-order expansion of Mc@Nlo 
and/or multi-jet merged runs are supported. The parton shower can be enabled in the run, but its 
dependency on the QCD input parameters can not be resolved yet as discussed in the previous 
section. The reweighting is then incomplete and the resulting prediction is not exact. Figure 1 
illustrates the validation for the internal reweighting approach. Shown is the PDF variation plot for 
the transverse-momentum distribution in Higgs-boson production at the LHC, comparing via a sin¬ 
gle run with reweighting, with dedicated event generation runs for the modified PDF sets. For this 
example, the reweighting run takes approximately 4 times longer than a usual run, but generates 
the same information (with less statistical fluctuation) as 52 dedicated runs. In this example (and 
all the following ones) the input PDF has been interfaced using LHAPDF 6 [11]. 
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Figure 1: Validation for the new internal reweighting feature of SHERPA 2.2.0. A PDF band for the CTIO 
set [12] is generated using a single run with reweighting enabled to generate all PDF error eigenvalues at the 
same time. This is compared with the band for the same PDF set, but now doing a dedicated event generation 
for each set, i.e. without any reweighting. 
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4. Automated Grid Creation using MCgrid 

The internal reweighting feature of Sherpa is useful for evaluating theory uncertainties for an 
observable. But when speed becomes very critical, because the number of re-evaluations needed 
becomes large (e.g. when fitting PDFs), more specialised methods must be used to reduce CPU time 
drastically. One established approach is based on cross section interpolation grids that represent a 
given differential cross section binned in the incoming partonic momentum fractions Xajb and the 
process scale jjL^ = jjLp = . There are two implementations for those grids called APPLgrid [4] 

and FAStNLO [5]. While interpolating grids are defined only for a specific final state observable, 
they allow for a very quick re-evaluation of the respective cross section and have considerably 
reduced disk space requirements, compared to full event data. Until recently, interfaces to fill these 
grids were only available for process-specialised code. This has changed, as two interfaces to 
multi-purpose Monte Carlo event generators are now available. 

One is called aMCfast [6]: It interfaces Lo and Nfo predictions using the FKS scheme [13] 
generated with MADGRAPH5_aMc@NLO [14] to APPLGRID. 

The second one is called MCGRID, which is implemented as a plugin to Rivet [10], such that 
(new or existing) Rivet analyses can be modified to create and fill interpolations grids. MCGRID 
reads in the events in the HepMC [9] event data format, such that in principle every Monte Carlo 
event generator can be used to produce grids via MCGRID. However, all required event weights 
must be inserted into the weight vector of the HepMC event. Currently, only the format that 
Sherpa writes into the weight vector is supported by MCgrid, which is based on the Catani- 
Seymour dipole subtraction scheme as illustrated in sec. 2. As for aMCfast, since the first release 
of MCgrid [7], APPLgrid interpolation grids can be filled with Lo and Nlo events. Here we 
report on two newly added features in the 2.0 release of MCGRID: the support of FAStNLO inter¬ 
polation grids and the exact treatment of the fixed-order expansion of Mc@Nlo events. This is the 
first time the FASTNLO package can be interfaced to a multi-purpose Monte Carlo event genera¬ 
tor. With APPLgrid and FAStNLO, all currently available interpolation tools for fixed-order QCD 
cross sections can now be used in conjunction with MCgrid. For more details on the revisions 
made in MCGRID 2.0, see [15]. Figure 2 illustrates for the case of inclusive jet production at LHC 
energies, that FASTNLO and APPLGRID interpolation grids created with MCGRID both reproduce 
the distribution of the observable with per-mille level accuracy, much smaller than experimental 
and inherent theoretical errors. 

5. Application: Residual PDF Dependencies of the Parton Shower 

As has been explored earlier [16], varying the QCD input parameters in parton showers can 
have a sizeable effect on the result. The possibility to correctly vary the full fixed-order expan¬ 
sion of an Mc@Nlo calculation—via the internal reweighting of Sherpa or using interpolation 
grids created with MCGRID —allows for the quantification of the residual dependencies of the par- 
ton showers that are beyond the fixed-order approximation. Focusing on the PDF being used, 
the method would go as follows: First, generate PDF error bands for a given observable for two 
different PDF sets A and B. Use the event generation for the central value member of A to ob¬ 
tain distributions for all members of B via internal reweighting or an interpolation grid. Now the 
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Figure 2: Validation for newly added support for FAStNLO grids in the 2.0 release of MCGRID. An inclusive 
jet Pt distribution is generated using a single run. At the same time, interpolation grids for FASTNLO and 
APPLgrid are written out via MCGRID. Their convolution with the same QCD input parameters used in the 
run are compared with the results from the direct run. The deviation is of a per-mille order. This has been 
achieved without any optimisations of the grid architecture (number of nodes/interpolation order). 

reweighted/interpolated PDF band can be compared to the directly generated band for B. Because 
the PDF has been varied everywhere but in the parton shower, any (statistically significant) devi¬ 
ation in this comparison between the two B bands would indicate the associated observable bin’s 
sensitivity to the PDF input of the parton shower. Finding sizeable deviations can motivate an im¬ 
plementation for tracking the dependencies of the parton shower. The sensitive bins would then 
serve as test cases to check if such an implementation works as expected, as the deviations must 
vanish then (or shrink if an approximation is used as e.g. only tracking the first parton shower emis¬ 
sion). Interpolation grids filled wifh addifional informalion for fhe parton shower would help in fhe 
filfing of PDFs fo dafa fhaf are nol appropriafely described by fixed-order calculations. 

In fig. 3 we exemplify such a sfudy wifh W boson production af fhe LHC, considering fhe pj 
disfribufion of fhe boson. The grid has been filled wifh residual dependencies on fhe PDF sef fhaf 
has been used during ifs creafion: fhe CTIO [12] cenfral value member. When convolufing if wifh 
fhe MSTW pdf sef [17] and comparing if wifh a sef of dedicafed runs, we observe deviafions in 
fhe low pr < 5 GeV region. This corresponds fo a shifl of fhe Sudakov peak posifion. The peak’s 
shape and posifion is sensifive fo mulfiple emissions and fo fhe behaviour of fhe PDFs af very small 
scales. Therefore a careful consideration of fhese sensifivifies musf precede a definife conclusion. 

6. Conclusion 

We have presented fwo novel mefhods for fasl re-evaluafions of fixed-order cross secfions 
suifable for fhe determination of fheorefical uncerfainfies or PDF/tti-fifs in an aufomafed way. Al- 
fhough bofh mefhods are based on reusing all fhe parfs of fhe calculafion fhaf are nol dependenl on 
QCD inpul parameters, Ihey do nof require evenl dafa fo be wrillen oul. 
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Figure 3: A Study of the residual PDF dependence of the parton shower of a grid that has been produced 
in an Mc@Nlo event generation with the central set of CTIO PDF. The MSTW-convolution of the grid 
deviates from the dedicated MSTW event generation in the low pr ^ 5 GeV region. 

The internal reweighting feature of Sherpa 2.2.0 is doing the reweighting on-the-fly, i.e. on 
an event-per-event basis. This is especially useful for simultaneously generating all distributions 
needed for specifying the theory errors for a given observable. MCGRID on the other hand is an in¬ 
terface for automated grid production using general-purpose event generators and Rivet. Because 
interpolation grids cut down the CPU time needed for a re-evaluation to the order of milliseconds, 
they are apt for PDF fitting. 

The support of both methods for filling the fixed-order expansion of an Mc@Nlo calculation 
gives rise to a way to study the dependencies of the parton shower on the QCD input parameters. 
We have exemplified this for the transverse momentum distribution of the W boson. 

In the future we plan to extend these studies, aiming to eventually implementing a method 
for taking into account (approximate) parton shower dependencies for the internal reweighting of 
Sherpa and for the creation of interpolation grids via MCgrid. Other features we want to add to 
MCgrid are Nnlo predictions and the support for multi-jet merging calculations. 
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